Critical Casimir force in ^He films: confirmation of finite-size scaling 
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We present new capacitance measurements of critical Casimir force-induced thinning of *He films 
near the superfluid/normal transition, focused on the region below T\ where the effect is the greatest. 
"He films of 238, 285, and 340 A thickness are adsorbed on N-doped silicon substrates with roughness 
« 8A. The Casimir force scaling function i9, deduced from the thinning of these three films, collapses 
onto a single universal curve, attaining a minimum i9 = — 1.30±0.03 at x — td^^'' — — 9.7±0.8A^''''. 
The collapse confirms the finite-size scaling origin of the dip in the film thickness. Separately, we 
also confirm the presence down to 2.13K of the Goldstone/surface fluctuation force, which makes 
the superfluid film '-^ 2A thinner than the normal film. 

PACS numbers: 68.35.Rh,64.60.Fr,67.40.Kh,67.70.+n 



An important focus in condensed matter physics is 
understanding how the properties of a thermodynamic 
system evolve as its size is shrunk to ever smaller 
dimensions. Near a continuous phase transition or 
critical point, the theory of finite-size scaling offers 
a testable prediction. According to finite-size scaling, 
the correction to the free energy per unit area of a pla- 
nar film of thickness d due to confinement of critical 
fiuctuations has a simple, universal form 
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where Tc is the transition temperature and the corre- 
lation length ^ = ^olil"'^ measures the spatial extent 
of fluctuations in the bulk, t = T/Tc — 1 is the reduced 
temperature. The scaling function 812 is predicted to 
be a dimensionless, universal fimction of the ratio d/^ 
and the boundary conditions that the order parameter 
satisfies at the confining interfaces. 

While finite-size scaling is applicable to all criti- 
cal systems, the most rigorous experimental tests to 
date have focused on the scaling behavior of the spe- 
cific heat anomaly of '*He near the superfluid transi- 
tion [^1^. This is due to the nearly- ideal, impurity- 
free nature of liquid '*He and the low-sensitivity of 
this system to gravitational rounding errors. For 
the superfluid transition, Tc = Tx = 2.17Q8K and 
ly = 0.67016 ± 0.00008 01 . For a 57 thick film, 
the magnitude and temperature dependence of the 
specific heat is found to be in reasonable agreement 
with finite-size scaling predictions Q. However, for 
films 500 — 7000A thick, the situation is not as clear- 
cut. The temperature-dependence of the specific heat 
is as expected from the universal d/^ dependence in 
Eq. Q. The maximum specific heat occurs at a com- 



mon value X = {d^o/S.y^" = td^'" = -9 ± iX^'"" for 
all films, where the negative x refers to the maximum 
occurring below T\. However, the magnitude of the 
specific heat shows an unexpected, systematic non- 
collapse (3]. 

The critical Casimir force is another fundamental 
manifestation of finite-size scaling that is open to ex- 
perimental testing. Just as the Casimir force be- 
tween two conducting plates arises due to the confine- 
ment of zero-point electromagnetic fiuctuations be- 
tween the plates , a completely analogous thermody- 
namic Casimir force is expected between the substrate 
and vapor interfaces of adsorbed liquid films, due 
to the confinement of critical fluctuations within the 
thickness of the film [111 110, [l^. The theoretically- 
predicted critical Casimir force per unit area 
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where Casimir scaling function {^{z) = 2Q{z) — 
zdQ/dz. Because in "^He films the superfiuid order 
parameter vanishes at both film interfaces, the criti- 
cal Casimir force is attractive (i? < 0) jsl, producing 
a dip in the equilibrium film thickness near Tx. The 
existence of this dip, first observed by 0, has been 
confirmed in a quantitative experiment using as sub- 
strates five pairs of capacitor plates made of polished 
Cu set at different heights above bulk liquid helium . 
The interpretation of this experiment is complicated 
by the roughness of the Cu surface, which changes the 
effective areas of the Cu plates and makes it impossi- 
ble to accurately determine the film thickness. AFM 
scans over 2500/im^ areas of the surfaces show they 
are not ideal, with 10-130 A rms roughness and oc- 
casional micron-deep scratches and dust particles, d 
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is calculated using DLP theory and by assuming 
that the Cu surfaces at different heights are flat. The 
result of the experimental analysis '8| is that the scal- 
ing function d exhibits a behavior suggestively similar 
to the specific heat. The minimum in the dip occurs 
at a common value x = — 9.2±0.2A^^'^ for all the films 
257-423 A thick. The temperature-dependence is ex- 
actly that expected from from the d/^ dependence in 
Eq. ©, but the magnitude of •& shows an unexpected 
non-collapse, the minimum of "d increasing systemat- 
ically from -1.85 to -1.4 as d increases from 257 A to 
423 A 15 . To address whether this systematic trend in 
the magnitude of i? is an artifact due to the non-ideal 
surface or is truly related to the non-collapse observed 
for the specific heat, we have undertaken improved ca- 
pacitance measurements of the critical Casimir force 
similar to but using flat N-doped silicon surfaces 
with roughness w 8 A. 
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FIG. 1: The sketch of the experimental cell. The level of 
the bulk *He at the bottom of the cell is measured using 
an annular capacitor (A). The *He film is adsorbed on 
the two silicon plates (B) attached to electrically-grounded 
copper guard rings (C) held 0.2 mm apart by Cu shim 
spacers. The temperature of the cell was measured using 
a germanium thermometer T attached to the top of the 
cell, which is calibrated vs. the ''He vapor curve and the 
lambda fixed point device (D) at the bottom. 

A sketch of the experimental cell machined from 
oxygen-free high conductivity Cu is presented in Fig. 
1. Two silicon (100) wafers highly-doped with phos- 
phorous (1-5 mV,/D ) are configured as parallel plates 
forming a capacitor with a gap G = 235/Ltm. Ac- 
cording to a tapping-mode AFM, the rms roughness 
of the wafers is 8 A. To minimize surface roughness, 
dust particles, and scratches, the experimental cell is 
washed, dried, assembled and sealed in the Penn Sate 
Nanofabrication facility, a class 10 clean room. Virgin 
wafers, completely intact, are rubber-cemented into 



the Cu guard rings used to position the electrodes. 
The top electrode is a 1 inch wafer, and the bottom 
is 2 inch. To minimize error from the fringe field, the 
top and bottom guard rings are grounded and the 1 
inch wafer is placed at virtual ground in the AC bridge 
circuit used to measure the capacitance C |^ ^| . To 
determine T\, a fixed point device anchored to the cell 
bottom is used, following the procedure described by 



The temperature control scheme of our experiment 
is similar to that of the original experiment of 0. 
A needle valve is used to close the helium fill line 
just above the cell. The data are taken with the 
cell slowly drifting through the lambda point, at 
10 — 40iiK/ h near T\ where equilibration takes longer 
and at 70 — 300 fiK/h below T\. We use two thermal 
control stages. The first outer stage is maintained 
at constant temperature with less than 50fj,K noise. 
To achieve a uniform temperature drift rate, we ap- 
ply heat to a second stage, just above the cell. After 
dosing helium into the cell, we typically observe signs 
of capillary condensation, where liquid droplets con- 
dense in the gap between the silicon electrodes. To get 
rid of these droplets, we very slowly (100 — 300fiK/h) 
thermally cycle the cell through T\, each time look- 
ing for a distinctive drop in capacitance that signals 
the flowing of liquid from the gap. This procedure 
is repeated until a reproducible C(T) dependence is 
obtained. 

To calculate the film thickness d from the measured 
C(T), we model C(T) as the equivalent capacitance 
due to three dielectric layers added in series: adsorbed 
film, vapor phase, and adsorbed film, obtaining 




where, if Co{T) is the temperature-dependent empty 
capacitance, the effective dielectric constant e{T) = 
C{T)/Co{T). As in 0, the dielectric constant of 
the film efiim = 1.05760 ± 0.00005 and the dielectric 
constant of the vapor tvapor is calculated using the 
Clausius-Mossotti equation, taking the molar polariz- 
ability of helium to be 0.123296 ± .OOOOiOcm^ / m.ol . 
The vapor density is calculated from the pressure 
P(T), using the second virial coefficient B{T) from 

m 

The temperature-dependence of Cq is due to a 
small linear increase in G caused by a combination 
of liquid surface tension acting on the Cu spacers 
and differential thermal contractions among the var- 
ious materials that make up the capacitor, includ- 
ing between the silicon wafer and the rubber cement 
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FIG. 2: The thinning of the film plotted vs. temperature 
T: (a) the data shown over a wide range of temperature, 
(b) a blow up near the minimum. The films are labeled 
by the thickness in A in the region above Tx where the 
Casimir force is negligible. 

underneath. In our data analysis, we assume that 
Co[T) = Co(rA)(l - 3.5 X 10-5(T - Tx)). This re- 
suits in a temperature-independent d for all films for 
T sufficiently above as well as below Tx- Each time 
we dose liquid into the cell to make a new film, we 
characteristically observe an additional small shift in 
G (and Co) on the order of 50ppm. To correct for this, 
we adjust Co(T\) in order to obtain the theoretically- 
predicted thickness in the regime above Tx where the 
critical Casimir force is negligible and the equilibrium 
thickness d on the silicon wafer is expected to be deter- 
mined solely by a competition between temperature- 
independent van der Waals and gravitational forces. 
In this regime, the film thickness is given by jj| 

where, on the left side mgh is the chemical potential 
due to gravity, fixed by the height h above the bulk 
liquid, where g is the gravitational acceleration and 
m the atomic weight of helium. On the right side 
is a simplified expression for the chemical potential 
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due to van der Waals forces, where 70 ~ 1950KA 
and di/2 ~ 2304 are substrate-specific interpolation 
parameters that characterize the net attraction of the 
helium to the silicon, including retardation effects ■ 
The parameters 70 and (ii/2 are approximate, ignoring 
the effect of the small 20 A natural oxide layer on the 
silicon. Nevertheless, the error is estimated to be less 
than 5% or 10 A and the same for all the films studied. 



T ' r 




FIG. 3: The scaling function vs the scaling variable x. 
The minimum occurs at a: = —9.7 ± O.SA^^'^, for all three 
films. 

In Fig. 121 we show the measured change in the film 
thickness in response to the temperature-dependent 
Casimir force near and below Tx, for three differ- 
ent values of the height h = 15.00, 8.01, and 4.22 
±0.05mm. The films are labeled by their thicknesses 
above Tx calculated from Eq. namely 238, 285, 
and 340 ±10A. As seen previously due to the 
Casimir force, thicker films exhibit larger dips which 
occur closer to Tx- Including the additional contribu- 
tion to the chemical potential from the critical Casimir 
force , the equilibrium film thickness is expected to 
be given by 

where V = 45.81 A /atom is the specific volume of 
liquid ^He, and d is the dimensionless scaling function 
for the Casimir force. The observed dip in d is due 
to -i? < 0, i.e. an attractive Casimir force between 
the substrate and vapor interfaces, as expected due 
to the superfluid order parameter satisfying Dirichlet 
boundary conditions at the two film interfaces p]- 

In Fig. 13 we show the Casimir force scaling function 
calculated using Eq. © and the data of Fig. [3 Be- 
cause it was necessary to disrupt data collection every 
2.5 days to transfer cryogen, each curve, which takes 
about 2 weeks to complete, actually consists of 4-5 
overlapping data sets that are spliced together; this 
results in additional noise and a small discrepancy 
very close to T\ not present in the earlier work [Sj. 
Nevertheless, within the scatter, all three data sets 
collapse onto a single curve, with a minimum value of 
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d = -1.30±0.03 at X = -gjiO.sA^^". The collapse 
of the data verifies that the dip in the film thickness 
near T\ is due to fluctuation-induced forces ||5|. ft 
is noteworthy that the measured shows quantitative 
agreement with the i5 obtained previously [^j for 423A 
thick ''He films on Cu, but disagrees with the results 
obtained for thinner films that, presumably, would be 
more sensitive to surface non-idealities. These results 
suggest the non-collapse is the result of inadequate 
corrections for the effects of surface roughness and not 
due to "i? depending on the additional off-coexistence 
variable hd^/" where A/v ^ 2.47 5^. This is expected 
to have important implications for the analysis of spe- 
cific heat and wetting experiments 0, 0, 0| . 

The new measurements, which focus on obtaining 
data near the minimum of the dip and over a wide 
range below Tx, confirm an additional important as- 
pect of earlier experiments: for all the films studied, 
we find the superfluid film is ~ 2A thinner than the 
normal film down to 2.13K. Experiments indicate that 
the onset of superfluidity in the films occurs some- 
where between x = -7 and -12A^^'^ |3,E3- Thus it 
has been suggested that the thinner superfluid film 
is caused by Casimir forces due to fluctuations in- 
volving superfluidity in the film, such as Goldstone 
modes, second and third sound 0|. As seen from 
Fig. 3, the thinning in the superfluid film is consis- 
tent with an asymptotic, low-temperature value of the 
Casimir force « -(0.30 ± 0.10)kBT/d^. This force 
is marginally larger than the — O.lSfcsr/d"^ force pre- 
dicted by [ig |. 

In summary, the current experiment confirms the 
validity of finite-size scaling formula for the critical 
Casimir force in adsorbed ^He films between 230 and 
340 A thick. Measurements down to ~ 2.13K also 
show the presence of an additional, non-critical, at- 
tractive fluctuation-induced force in the superfluid 
fllm. Our study underscores the importance of smooth 
surfaces for these types of measurements. For fu- 
ture work, it would be desirable to test the scaling of 
Casimir forces in a much wider thickness range that 
overlaps the range covered by specific heat measure- 
ments. We would like to thank A. Maciolek, R. Zandi, 
J. Rudnick, M. Krech, D. Dantchev, S. Dietrich, M. 
Kardar, S. Balibar, G. Williams, F. M. Gasparini, and 
M. W. Cole for useful discussions and comments. This 
research was supported by NASA grant NAG3-2891. 
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